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ABSTRACT 


There existed aneed for an interactive program hat 
would provide the usar assistant? in solving applications of 
damnear control thesry. The Liasar Control Program (LINCON) 
and its user's guide satisfy this need. A series of ten 
interactive programs are presented which permit the user to 
carry out analysis, design and sSimulatioa of a broad class 
Sepeenear ccitrol problems. 
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transfer function aad tine response programs; and modern 


controls programs. Examples for aach are work2d within each 
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terminal session section. 
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I. INTRODUCTION 


The purpose of this thesis was to update an existing 
program which provid2s assistance in solving computational 
problems associated with the stady and application of linear 


COncrol theory. tiewkanear Cone rols Program ({LINCON) was 


th 


irst developed by Msisa f 1] and adapted for batch use at 
NPS by Desjardins [ 2]. Miata ne wast ginal wnten* of 
this thesis was simply to take Desjardins' adapted version 
of Melsa's LINCON and further adapt it by making it 
interactive, LINCON soon began to grow as other routines 
were incorporated, 4s will be noted, until its present forn 
was achieved. 

LINCON, as such, is @ higa level applications software 
sysctem made up of 32 large number orf program tools for 


interactive analysis, design ani simulation of a broad class 


Mee linear control problems. Weta EL Ne ON, users can 
Bemcen=erate on th2ic specializ3d applications rather than 
System design and routine program development, thereby 


Saving valuable tins. 


ies essumedq that the reader is familiar with the b 


pv 
"i 
} 
©) 


it 
i4 
O 
= | 


BeeceptS 929i iinear control theory as may pe obtained 


it, 


any one o @ numb2r> of available <textbooks (see the 


ee) 





biblioagraph y). WS suena, thea LINCON subprograms are 
presented in a user- sclented fashion. First, their purpose 
and some general rules that may apply are given; then the 
input requirements are present2i and the expected outputs 
are described. Examples for sath are work2d cut including a 


copy of the terminal session and the final results. 
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II. SYSTEM IVERVIEW 


During LINCON 's latest development, underlying 


guidelines called for concepts which accounted for the needs 


of the computer syst2m, the programmer, and the user. 


The guidelines followed jiuring this latest development 


stage were: 


(1) 


(2) 


(3) 


(4) 


Operation of the system should be in an on-line 
interactive node such that data can be easily input 
to (or output from) ths system and readily acc2ssed 
POtOVenbatlecatLOn, €Xamination, and processing. 


Program development should be in a high-level 


th 


language in order ae. acilitate software 
implementatioa and promote machine independence. 
The software should b2 nodular in structurs so chat 


programs 7ca2 b2 modified or TPSerled wwe = Nout 


(Dp 


@ebece ing @xisSting programs. 
Programs should be iavoked by néans of logical 


Taras] = on 


J 


procedures or ccommanis which mininize i 


time and whith are user-oriented s9 that people car 


} te 
it 
He Ie 
VU 
Ee 
vf 
th 
|- 
14 
" 
ct 
OE 
MD 
QD 
© 
= 
|. 
3 
iQ 
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O 
=| 
“fC 
jo 
ck 
(Dd 
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Operate the system w: 


experts. 


Wy 





fe «LNT ERACTIVE OPERATION 
PRremeepolsecacsons OC. LENCIN are focused on interactive 
processing. Experiznce has shown that interactive on-line 


communication has many advantag2sS in a research environment 


because it offers the opportunity to make observations and 
select alternate courses of action in a more flexible manner 
than with batch processing { 3}. LINCON is organized around 
a collection cf inter-relatad command programs, each of 
which performs a sp2tified funstion and can be executed by 


means of a simple keyboard initiation sequénce. 


Be HIGH LEVEL LANSUASE 
Mimemporcanrt Leature in the d@¢sign of LINCON is thak it 


was inpl 


(D 
=) 


ented 19 a high level language. Program 
Gevelopment in assemolLy languag2 is mors time consuming and 
results in system dependent software. 

LINCON is programmed in ANSI standard FORTRAN and 
BebtOows the conventions of FORTRAN IV. FORTRAN has been 


found t9 be a useful Language for several reasons: 


(1) Since some form of FOXTRAN is ivailable on most 
GOMPUTETS, DENeet is haghiy poetabie from one 
computer to another. A FORTRAN based system is 


nel prul Leb importing progtams as well as expor*in 
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them. Of course, FORTRAN compilers don't all follow 


the Same standards so there Gan. “Stall be 
dest ficulties. 

(2) FORTRAN is a simple 23n ough language that relatively 
complex programs can b2 implemanted in a short 
period of time. Most scientific and research 
personnel know FIRTRAN sufficiently well to write 
their own programs if n2cessary. 

(3) Algorithms can be test2i and implemented in FORTRAN 


and later convertsd to assembly language versions if 


more speed and =fficiincy are nscessary. ayia 
procedure has the farther Demertit of aiding 
portability such that even if parts have been 


converted to assembler, equivalent FORTRAN versions 


are availapbls. 


e. MODULAR SOrTWARE 

Overall system flexibilty is achieved by means of 
modularity. LINCON is actually nade up of 2 large number of 
independent command programs. sach ccmmand program stands 


ees )«©6COW ET hUWWAith Uhthe «ability 3 t2 take sont form of inodut, 


possibly supplied by some oprsvisus command, and generate 
SCme Crorm or output, possibly *o be us2d by a follow-up 
ccmmand. 
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D. USER-ORIENTED JPER ATION : 

An important aspact in this modification of LINCON was 
to make the commands uSér-orianted so that operating the 
system does not r2quire an engineering or computer 
background. This was achisv2i with a standard terminal 
keyboard by using phonetic characters which relate ‘*9 the 
function which the commani is to perform. A good 
Semoinetion of brevity and clarity is built into LINCON to 
avoid having to push 2xtra buttons on the keyboard. whils act 


the same time preventing ambiguity. 


Bee oNPUT RESTRICTIONS AND LIMITATIONS 


Although the input requirenents are fully described in 
the presentation of 2ach progran, there are several input 
format similarities ased by all of them. For ease of use, 


and, more honestly, for ¢ase of programming, most of the 
daza input is groupei in the san2 arrangement. 

meee rlmse iWput 26 every program is used to identify the 
Meemetem for references and for dvutput dat2. A naxinum of 
twenty alphaenumeric characters can be used. Tyas 
Testriction was not 32 system limitation but a programmer 


decision. 





The next input common to all programs is the dimension 
Seeene plant matrix >t A matrix. The format is 11 which 
would normally restrist the user to 2 maximum matrix sizs¢ of 
9x9, however, du2 to, again, a programmer decision, a 
dimension size not to exceed 8 is requested. The reasoning 
behind this was due, in part, tor the printer. The NPS 
printer is capable of printiny 133 characters on a line. 
Since the output format to the printer 1s 8516.6 this 
naturally limits one to @ numbers per line. Six places are 
rormally considerei ascessary f5r good accuracy. A solution 
that would have iead to an unlinited matrix size would have 
Beer to incorporate a "wrap around" routine within the 
program. After attempting this it was decided the results 


Bepeer yuSt <0oD difficult to read. 


fu 


Matrices are ent2red on 


(I) 


alement az a tim? beginning 
Meeer Clement 1,1 and Sontinuing across the row. The next row 
is then entered, and the process continues until all 


elements have been entered. After the matrix is entered, the 


SempletS Matrix is automatically brought to the screan for 


Ho 


review and possible correction. If a chanya to the matrix is 
desired the user singly enters the row number and column 


number without a separating conna. For example, 35 wouléd 





indicate the element in row 3 and column 5. REwersbsing 
prompted the change? is entered. A review of the matrix is 
again screened. The aser 1S again prompt2d for any possible 
changes. This procedure continues until all changes have 
been made. 

Any special requirements or limitations will be brought 


to the user's attention within 2ach prograa presentation. 
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Ae INTRODUCTION 
In this chapter four prograns are discussed which may be 
used for the analysis and design of linear control systems 


represented in state variable form as 


x(t) = Ax(t) + bu(t) (3. A-1) 
w(t) = KE r(t) - k x(t) ] (3. A-2) 
a 
y(t) = ¢€ X(t) (3. A-3) 
mone first, Basmat, as “ta= Basic Matrix manipulation 
program which is ased Co= comping the decerminaat, 


inverse, characteristic polynonial, and @igéenvaiues for 4a 
Square matrix A. Eyeaddisce on, BASMAT@WG1 1 calculate the 
mare transition matrix and the PHI(s) matrix. The second 
Meegram, PRFEXP, calculates the partial fraction expansion 
Sea polynomial. Piew chi ed pEotJtam, ROOTS, calculates the 
Meets Of 2 polynomial. Mac pecan program 2S used. £95 
MeeecMining the tims respons? of linear control systems. 


RTRESP will determins the ratis1al time response of a system 


ct 


Mee. Osed-Lorm proviied that the input function r(t) has a 


ya 





rational time response and that there be no repeated 
eigenvalues in the combination 3f the system and input. re 
should be noted her2, and will be again in the actual 
discussion cf the program, that by setting c(t) and K equal 
mo Zero, unforced 22d open-l2op systems may be studied, 


respectfully. 
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Be BASIC MATRIX PROGRAM ( SASMAT) 


the 
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lan) 


Given plant matrix A, BASMAT an 


men .OWLNnGg: 


(1) the determinant of A, dat A 


(2) inverse of A, A-~? 


(3) Skhezactere seis polyaomial, det(slI-A) 


(4) eigenvalues of th2 tharacteristic po 


State transition matrix, ® (+) 


>) 


PHI(s) matrix, ®(s) 


(6) (sI-A) 7? 
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eSaSMAT}). 
ANALYSIS IS COMPLETE. DO YOU WANT TO RUN LINCON AGAIN? 
1g 


LINCON IS NOW TERY INA ED. 


Pies eSOMDesm Lesulits are shown in Apoendix A. 
Meee pretation of the detemmin2rt, inverse and eigenvalues 
Meee ScraAightrorward. The PHI(s) Nese= eso a COMDLNStL On oF 


meme numerator coefficients ani characteristic polynomial. 


The results can be interpreted is 


Page 





1 "s2¢5-1.5 0.0 0 
$ (5) = se eee enn ---- F s2-2.0st1 s-1 
s2-2.585+1.5 ; ie, ~0.5s+0.5 S?+s~2 
fed, Simiarily, the fifth term of the state transition 
matrix is 
® (+) = -0.06/7exp().82t) +1.0/72xp (-1.82t) -2.26x10-Sexp (t) 


ZS 





~_ 


' | a 
Pe abst . Oe 
#203, bore te Es ove | 


4 


aac . 


te o3* 30 WSS qe oESnee 


fo, @429¢ ygeshatueee 





Se PARTIAL FRACTION EXPANSION PROGRAM (PRFEXP) 
PRFEXP can calculate a partial fraction expansion given 


mecaerOndl ratio of two polynomliais in the forn 
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DO YOU WANT TO RUN THE PROGRAM AGAIN? 


Muaeroks LS COMPLETE. DO™YOU WANT TO RUN LENCON AGAIN? 


LINCON IS NOW TERMINATED. 


The results are shown in Appeniix 8B. Tieee=nDretatzon OF 


these results are: 
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De POLYNOMIAL ROOTS PROGRAM (RIITS) 


ROOTS can calculate the roots of a polynomiai of degree 


less than or equal t> eight. Siv2n a pelynomial in the form 


Pegs) = "py * DLS + pyse + way + s’ (3. D-1) 


the coefficient of s* must be unity. fhe program ouput 


lists the polynomial coefficients for reference and <he real 
and imaginary roots. 
1. Terminal Session Exampls 


The following polynomial is to be factored: 


P(s) = -5s3 + 15s2 - 25.585 + JU 


Putting the polynomi2l in the required forn yielis 


P(s) = s3 - 352 + 5.15 - 0.8 
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aeons Uo TO FIND THE RIOTS OF A POLYNOMIAL 
DEGREES S SHHAN JR BQUAL TO EIGHT. 
Pies  eNrene tas PROBLEM. TDENPITFICALTION 
(NOT TO EXCEED 2) CHARACTERS). 
thesis exarple 
3 Peet as ORDER JP THE SYSTE2M (UP TO 8). 
ROOTS OF A POLYNOMIAL 
Meron LHe PORYNOALAL COEFFICIENTS IN ASCENDING 
ORDER OF 5S. 
PEMA ENG==-THE HISHEST ORDER COEFFICIENT MUST BE 
COEFF(1)= 
26 
COEFF(2)= 
Bis 1 
COEFF(3)= 
5. 
, COEFF(4)= 
Poti Nomen becOnrr LT ClENTS - LN ASCENDINS POWERS O 
Eevee -Ut oO. TUETO0 “3S. J0ETO0  T.00E+00 
THE ROOTS ARE REAL PART IMAG. PART 
=-T4iet00 1.62E+0) 
~1.41E+00 -1.62E+#+00 
lesbo ol O.0 
Mae SoeCONGCLUDES THE POLYNOMIAL ROOTS PROGRAM. 
DO YOU WANT TO RUN THE PROGRAM AGAIN? 
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minions £5 COMPLETES DO YOU WANT TO RUN LINCON 


LINCON IS NOW _ TERMINATED. | 
The results shown 12 Appendix C are 


interpreted. 
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Be RATIONAL TIME RESPONSE PROGRAM (RTRESP) 

The time respons2 in clos2d form cf a linear control 
system described by Egs. (3.A-1), (3.A7-2) and (3.A~-3) is 
calculated by this oprogranm. The user must define the 


initial conditions x()) and the rational Laplace transform 


ct 


Bumeeye Scalar forcing Eunction rit). The theoretical 
concepts involved in the development of the computer codes 
are described by Melsa { 1]. 


1. Input Restrictions And Limitatiors 
Preece ene E£OLLOWing in & format (floating point) : 
1) the elements of the plant matrix A 
Bie eno CONT TOL Vostor 'b 
Bye ene OUtpUs Vector < 
4) «he feed backback co2afficent vector k 
Sieene comerolLler gain < 


and 


Seecie ©6pt22i “Conditions vector x(0). 
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so that the coeffici2ants of s and s’ are each unity with 
mee g 2 0. Due to programming Iimitations it is necessary 
that the order of th2 system plus the dimension of D(s) be 
less than or equal t9 ten. 
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memthe Laplace transform of the input function is 


a (6a 


Smamcne initial conditions ars zero. 

Note here that the open-loop rational time response 
mey be calculated by setting the feedback coefficient vector 
Mmeeguel ¢0 a zero vector. Also, SG (Orly ah 2naekes. 
Semel czOn response is desired, the input function gain G is 


Bert tO zero. 


fence on 
SxecULION BEGINS... : 
LINCON CONSISTS OF THE FOLLOWINGS SUBPROGRAMS: 
BASIC MATRIX MANIPULATION - <BASMAT> 
RATIONAL TIME RESPONSE - <RIRESP> 
STATE VARIABLE FEEDBACK - <STVAR> 
CONTROLLABILITY AND OBSERVA3ILITY - <OBSCON> 
LUENBERGER OBSERVER - <LUEND> 
OPTIMAL CONTROL/KALMAN FILTERS - <RICAIID 
DISCRETE TIME KALMAN FILTER - <KALMAND> 
OPTIMAL CONTROL - <OPTCOND> 
PARTIAL FRACTION EXPANSION - <PRFEXP> 
ROOTS OF A POLYNOMIAL - <RODTS> 
Memeo ONE OF THE SUBPROGAAMS ENTER THE NAME BETWEEN THE 
SYMBOLS < >. 
riresp 
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EN@TER De EEMENTS OF THE JUTPOT VECTOR--C. 
C(1)= 
0 ( 
C(2)= 
4 (2) 
C(3)= 
0 
Mite oat R nx (OUTPUT VECTOR} 
1. COE+00 
GY“. QOE-01 
1. 0O0E+00 
DOmYOU ~GANT TO CHANGE ANY BELEMENT(S) OF THE MATRIX? 
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DO £0U WANT TO CHANGE THE VALUE OF THE SAIN? 


RO}. 2 eee OF “THE ZNITLAL CONDITIONS VECTOR -- 
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0.0 a 
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merce at CONDITIONS VECTOR-X (0) 
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DO YOU WANT TO CHANGE THE VALUE OF THE GAIN? 


ENTER THE NUMERATOR BY POLYNOMIAL 
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FACTORED ROOT FORM. 
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peteoroiS TS COMPLETE. 


Ee NC 


ESP) 


ON, 25 »NOW TERMI 


PROGRAM 


DOR Gmewatt tO RUN LINCON ASALN? 


Nat cD. 





SR ODeEnG:x PD are 


shown 


results 


computer 


The 


interpreted as 


0.0 


x, (t) 


VeseeneGiea a=) 


be 
a 


p(-0.5 


a7 





iV. 


= 


ODERN CONTROL PROGRAMS 


A. INTRODUCTION 

In this chapter six programs are presented and discussed 
which may be used for the analysis and design of control 
systems. 

meeot. rst, OBSCIN, is used to find the observability 
index and ccntrollability of a system. The next two programs 
are used +o0design <alman filters. RICAII and KALHAN can 
merant he feedback ani controi jains necessary to optimize a 
Benetton for €lither continudus or discrete systems, 
respectively. The iast three programs may be used to disign 
@eewemal iz=near control svstems. SIVAR is particularly us?2ful 
Mameecne design of liarear stat2 variable feedback control 
sys ems. Zi o6may ~«©6p2)6hlUUSed t>06Saiculacte both open- and 


Meesea-LOop <tranSt2r functions and also nas the ability <o 


Meseqn a closed-ls0p system Erom desired transfer function 
Beecartications. LUEN LS used to design a “Gomb=ined, 
reduced-crder, observer-controlisr to achieve a desired 
Crcsed-loop transfer function from a system where some of 
THE States are inaccessible. OPTCON will minimize a given 


ese c£uncticn producing @ scalac control. 
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Dee OBSERVABILTIY INDEX AND CONTROLLABILITY PROGRAM 


This 


madex and centrollability of th: 


x(t) = Ax (t) 
y(t) = Cx(t 


where 


a 
tl 


ic 
i 


program is 


State vectoc 


t> determine the 


linear systen 


(n-dimensional vector) 


centrol vector (l-dimensional vector) 


Ve OuLpw=] VSCtOr (N=d2menasional vector) 


Feo 
T 
re) 
ra 
fa 
= 


joo 
tl 
3 


C= aX 2 Maer x 


The sbservability 


integer such that ths 


T T 
CC, A Crees, (A JT 1C) 


mos rank n. 


Yomia. Mat. > iat 


roel) < 2 


The above systen 


Siew esen POSS=aDLe, 


BPS, 


(OBS CON) 


observability 


(4. B-1) 


(4. B-2) 


momsaed 26 D= “CORtzoOllabpi 


ho mininun 





Mmcenstrained control vector, to Lorce the system from ar 
initial state of x(0) to some other stat2 in a finite time 
interval [ 4]. 


Tf the user desires just the observability index to be 


Sereulated, enter tha B matrix 1s a zero matrix. Likewise, 
if just the controllability of tha system is desired, enter 
mee C Matrix as 2 z6ér> matrix. The unobservaple or 


UW 
(D 
— 
Va 
qc 
aie 
iD 
ee | 
O 
rh 


uncontrollable system respon SOurse Lgnered. 
Ws Input Restrict ions 
Other than the limitations cS the problen 
m@emere2cati0n, (20 characters) and system size (8x8) thers 
Meee NO restrictions +o OBSCON. A reminder here may  be- 


helpful, however. Remember £5 snter your elements of the 


Weeest ces in F format. 


~ 


Pomeloe foeat S25520n Exanpl= 
Mies section COnteias 2 Session for a specific 
oreblen. Commands 2atered by th2 user are in lower case. 


(D 
cr 
( 
(i) 
ct 
‘D 
Qu 


The following system is *95 b 


MORO) 1.0 9 0.6 0 aaa 
mee y = -1.0 -0.5 0 Bete) OO 10 a (+) 
920 029 io 050 0.0 ° 
p | 


+9 





~ncon 


e 
2 


1 


BETWEEN THE 


“~ 

ee = fr] 
2) O = 
0 O A et 
foo mm TH a 
rr Mm -—FH 
O O ata fx 
O V OF A 
mA Fetuoof Ft 
Aue4 © t Mma PSs 
MO «<< VA fe 
DEAA> GQ fu fa 
MMGe Ih wm H 

GU) eGHI VY AA 
OMS WH Vw 
ZV/UGRH mI FH 
Hi RHinmaA ts ean 
wl maVeaaeia Ow 
QO Vi Sem ama wais 


12M 53 Aue 
MmOAAMAGYV FIM 
MUA esheets! fst 
$4 UU) F2] crt Be 1) 1 Qin 
efits MND we 


Z=HNSOQAQHRAUN, 
HHH H BH 
OHMmMAHMMOHOMAO 
IW ghU O fu 
(Lt at Gs) 


Nadal INDEX AND 
iC A EON 


G 
i! 
PLA? EYAL RIX =A". 


+1 


Poo or Las 


fe} fr,ft] « 


OMAP AMABHHAHOBA 


eal ae. 
>< O= 
fx) E+ pb 

— 


Pebicke lee ORDER OF THE sY¥STcam (UP TO 8). 


WNTER THE 
A (1,1) 
Cty, 2) 
ABC ¢ 3) 
A(2,1)= 

A (2,3) 
Bs » 1) 
A(3,2)= 
OS, 3) 


NOW, 


® 





Pas A MATR LX NGO MATRIX) 
1.00E#*00 0.0 
-~1.00E+00 -5. QO0E-01 1.00E+09 
: 0.0 1. OOE+00 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
BYloR THE NOMBER OF COLUNNS OF THE B MATRIX. 
stn 
2.0 
a (1,2) = 


= BLSMENTS OF THE B MATRIX. 


As 


0 ae 


B (2,2) = 
ose 4) 


B(3,1)= 
oo or" 


B (3,2) = 
A” | ) 


THE B MATRIX 


2-Q0E+00 1.QQ0E+00 
0.0 1. QOE+00 


DO YOU WANT "TO CHANGE ANY ELEMENT OF THE MATRIX? 
ieee IE NOABER JF OUTPUIS. 


THE SLEMENTS OF The C WATRIX. 


+2 





THE C MATRIX 
1.00E+00 0.0 

OE+00 1.00E+00 0.0 

OE+00 1.00E+00 0.90 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


OBSERVABILITY INDEX = 2 


THE SYSTEM (A, 3) IS UNCONTRILLABLE 


tok ook ok ak ok tok Zc kk ako ak kkk ok go kk tak kk kak kk ae doko kek ok ok ak 2k 
DO YOU WANT TO RJN OBSCON ASAIN? 
n 


THIS CONCLUDES THE OBSERVABILITY INDEX/CONTROLLABILI TY 
PROGRAM (OBSCON) 
ANALYSIS IS COMPLETE. DO YOJ WANT TO RUN LINCON AGAIN? 


LINCON IS NOW [TERMINAY ED. 
The complete results pr2sented in Appendix E should 


be self-explanatory t> the user. 
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fe OPTIMAL CONTROL/KALMAN FILTER PROGRAM (RICATTI) 

The transient solution, matrix gains, to the Riccat: 
differential egquatioas for th2 state-regulator controller 
and the continuous Kalman filter are calculated by the 
RICATI program. 

Given a state-ragulator problem with the linear,tine- 


invariant system 


x(t) = Ax(t) + Bu (+) (4..C-1) 
men = Cx (*) (4. C=2) 


mien kiccat: equation is defined as 
——_ = 
Peeepee- “E(C)A=A P(t) +2 (f) BR-*B 2 (t) -Q (4. C-35) 


ween P (to) aS the DOURndery coad2ztion. Pe eu (ie) is not 


@emetrained, a gain natrix can be found such that the cost 


mune t On 
ty 
J=1/2( x" (tC, JB(Te) X(t) ev2[tets Qx (t) tu’ (t) Ru(-) jade (4. C-4) 
to 
mememenemized { Sj. Such 3 gain natrix is isfined as 
G.(t) = B-1B'B (t) (4. 0-5) 
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w(t) = -S, (t)x (t) (4. C-6) 


fre. transient solutidn is solvai by the computer. Notes that 
the output from the somputer for the gain matrix does not 
inciude the negative sign of th> feedback loop. 

Given a continous Kalman filter problen with the linear, 


time-invariant systen 


x (t) 


Tt) 
Tce 
j>< 
am 

ct 

a 


Bw (t) (4. C-7) 


iN 
‘Ta 
~~ 
EF 
ap 


y¥€&) = v(t) (4. C-8) 


where w(t), the random process forcing input and v(t), ‘the 
measurement noise, h2ve covariaice matrices of R, the random 
mpu: CcOVarliance matrix, ani 9, the neasurement noise 
covariance matrix, HeSpecuPuly,nstne Riccaci equation is 


defined as 


P(t) = AB(t)+B(+)A +BQB -B(t)C 2-!CB (+) (4. C-9) 
wah 
Blt) = ECL R(t.) -¥ ko) -E 1 (ie =10) 


Gemeene 210it2al condition boundasy. The g2in matrix found is 


aefined as 


t) = R-tcp(t) Gao 





1. Pepwinal Session 


example presented by Melsa [{ 1]. 


Exampl 


+ion contains 


iw 


an SESS OV LOnreweme — SpeCiiic 


Commands entered by the 


user are in lower cas@. Given tha second order linear system 


e t) 1 1.0 ty 
xX = x 
x ( 0.5 0.0 =! 
| 
: eo 050° 7 
= zie 
x(t) 0.0 2.0 =! 


determine the optimal transient 


gains for 


cr s| 
_ ile ie 
c= 2 = 
Vs 9, 
& j 
fr | 
= 0.0 
R=R= 
0.0 5 
L j 
Mes the control option an 
meme Or 10.0 is used. For 
Samu.) and 4 final tine of 
eeeon, the initial condi 


Non Mace sx 


1.0 9.0 
. u(t) 
Van 140, 


Eespense control and filter 
Mimic amenor O50 and 3. final 
Biewe ee OD rion ah Initial tine 
Mees eUSeads Also; LOz “the fi leer 


TS) chosen £0 be 





on 


and 


ten equally soaced values of G 


an~enval t 


mor DOth Options 
the time 
eaval to 10. 
Pamicon 
PAeCUTIO 


Tt 
pw 


ort 
E+ 
>, 
fx 
{x 
4 
aN 4 
ee i, mM 
Ww) © 
= Sy 25 fry 
Pi WW HH = 
mK m Ft <x 
8) OQ tA a} 
O Vv UB A 
lA Hee a fd 
Qu E4 t (os 4 oe 
M et VH ww B&B 
MZAA D> eq fry 
WWMM jt Mm & 
co ett Vv GAR 
INNS WH Wwe 
m\/ FI mt Fi, 
Mt FHINMAG) § fr 
Twi mV sO Se 


HH I MmHneZove 
HOt MYawHNAOHWH «ae 
OH MMV HUW 1% 
S.F43()M wae, 09 
TWO’ at UeHo 
ho=emM sO 
MmOANMGMAAKV HA 
HORS hleagst fIsim 
Ht WF] crt & he oS 
ofa Aly MN BmaWN 
COTM Nhs aIOm 


OF tre YAMA HYHOnM 


20 2% ONGaRtert AV 


VORMmMS GaSHMO 
HHHEISHOHAH 
WH 

eetHODMHUtOno,) 
MMNUHNOACAaMmDmA 


O 
fH 


LINCON 


Po Boa RE RI ASA. 
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aoe CONTROL GNPUTS. 


(PLANT 4ATRIX) 


—2. VOE+00 


Bienoni Ss) OF 


MATRIX 

0.9 

WANT TO CHANGE ANY 
THE NUMB 


Penne ORDER OF THE SYSTEM (UP TO 8). 


cZPI BMH MORTAR I 


NTER THE 
NTER 


DO YOU 


NOW, 
A (2,1) 


£2iCc 





Pet is EL sNro OF TLHE DESTREBUTION MATRIX--B. 


B(1,1) = 
10 UO, 


1,2)= 
mo" 


B(2,1)= 
0.0 


B(2,2)= 
1.0 


mae BS ONAL RIX (DISTRISUTION MATRIX) 


imeoors 00 0.0 
0.0 1. OOE+00 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
aiwon THE NUMBER 9F OBSERVASLE OUTPUTS. 


Peon LHe ELEMENTS OF THE MEASUREMENT MATRIX--C. 


ert) = 
0 (1, 1) 


G(1,.2) = 
0 (1, 2) 


C(2,1)= 
0.0 OA 


C(2,2)= 
0 (2,2) 


ms CMADREX (MEASUREMENT 4ATRI<Z) 


ie COE+OOm 9.9 
0.0 2. 00E+00 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


h 
YOU HAVE TWO OPTIINS AVAILABLE: 
(1) THE CONTROL OPTION FIR SOLVING STATE-REGULATOR 
PROBLEMS, IR 
(2) THE FILTER OPTION FOR SOLVING A CONTINUOUS KALMAN 
FILTER PROBLEM. 
FOR THE CONTROL DJPTION, ENTER THE LETTER C. 
FOR THE FILTER OPTION, ENTER THE LETTER F. 
@e 
ENTER THE ELEMENTS OF THE CONTROL WEIGHTING MATRIX--2. 
R(i, 1) = 
1.0 iG 
4) 
R(2,1)= 
i |OCO 
Ri2,2)= 
>. iad 
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Paden GBATRIX (CONTXOL WEIGATING MATRIX) 


1.90E+00 0.0 
0.0 2. OOE+00 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


tipo LAS ELEMENTS OF THE SPATE WEIGHTING MATRIX~~Q. 


1,1) = 
1.0"! ) 


1,2)= 
a 


2,1)= 
1.0°§ ) 


2,2) = 
1.02! ) 


mae e MATRGX (STATE WELGHITNG MATRIX) 


1,.Q0Q0E+00 1. QQE+900 
1.00E+00 1. 00E+00 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


J 


<a Mitte AL Pee esOn THE TRANSIENT RESPONSE. 


10 saa Meee TEME FOR THS TRANSZENT RESPONSE. 
; ee Peete oman oleND RESSONS = 
10 


WE ie EE eK ie ae ie Re KK oko OK Kk ok Kk 


meow CONTROL OPTION *** 
feok kok ok ok tok ok ook eek ok Kk ok ok kk 


ENTER THE SLEMENTS OF THE TERMINAL BOUNDARY VALUE 


MATRIX-->. 
P4i,.t)= 
oR (1, 1) 
D(1,2)= 
a2 0 Se 
2(2,1)= 
0.0 (2, t) 
P(2,2)= 
0.0 (2, 2) 
THE P MATRIX (TERMINAL BOUNDARY VALUE MATRIX) 
Oe 0 0.0 
0:0 0.0 


Dem OU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
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TRANSIENT SOLUTION 


1. OOOE+01 


9. OOOE+00 
E 
E 


Oo Orc 
. ©€@ e 6? 
, a aD Nc 
tnnre 

vaDOQ 

Foti f 

So «Gfx 
HUOWS 

Oe He "st 
e @ e #$¢ 
Coe Oc 


8. 000E+00 
E 
i 


7. Q00E+0) 
E 
E 


6. Q00E+0) 
E 
E 


5. 000E+0) 
a 
ie 


4. Q000Et+0) 


3. OOOE+ 090 


2. 000Et+0) 


1. OOOE+GO) 


2G eNo- U6 


We KK ee ee KE 2 Ke oe oe ee ee KK KO eK KO OK OK KOK KK OK KKK OK KK 





DOMVOUMWANT THE FILTER OPLIIN? 

Teens FomeGnanGeo TO THE A, B, OR © MATRICES? 

Rieepe lH PRE BEVCHANGES TO THE R OR QO MATRICES? 
ee THE INITIAL TIME FOR [THE TRANSIENT RESPONSE. 
eo Tage Stab PME FOR THE “TRANSZIENT RESPONSE. 


THe WUMeeR OF POINTS JF THE TRANSIENT RESPONSE TO 
“i N (<109) 


yok gok ko doko kk ok fk 4 kok 


Powe. Usb LON | eae 
ek ek ok ok tok aa OK KK KK 


Peete LHe ELEMENTS OF THE INITIAL BOUNDARY VALUE 
Balaki~- -P. 


ei) = 
ao : 


P(1,2)= 
noe (1,2) 


P(2,1)= 
0.0 (2, 1) 


THE P MATRIX (INITIAL BOJNDARY VALUE MATRIX) 


0.0 0.0 
0.0 0.0 


Peer oaus WANT TO CHANGE ANY ELEMENT Or THE MATRIX? 


TRANSIENT SOLUTION 


TIME = Q.0 
GAINS 
0.9 0.9 
QO. O58. 


NO 
+3 
twaH 
!'’rHo 
OOo-24 


@ 
e 8 
tzid ire tx 


OW 
——/) j| 


COND 

@ 6 
ins 
kas GD et 
Cita 
1 ttt 
OO's 
—_a—(/) ji 


tJ 


ho 
e686 
n~] 





TIME = - QO0E+00 
GAINS 
2. /8E-01 2e13E-91 
8.S3E-01 7.06E-01 
TIME = 2.500E+00 
GAINS 
2. /79E~-01 2- 13E-01 
8.53E-01 7.06E-01 
TIME = - OO0E+00 
GAINS 
2.79F-01 2.13E-01 
8.53E~01 7. 06E-01 
TIME = 3.500E+0090 
GAINS 
2279E-01 2.72E-01 
8.53E-01 7.06E-01 
TIME = 4&.O00E+09 
GAINS 
2e-79E-01 2.13E-01 
8.53E-01 7.06E-01 
TIME = 4.500E+00 
GAINS 
2-79E~01 2.12E-01 
8.53E-01 7. 062-01 
TIME = 5.0Q000E+09 
GALNS 
2e79E~01 2.132-01 
8.53E~-01 7.06E-01 


He ake 0 ok keke keke kk ok fk ke ak ok ke ok ok ak kk ok ke ok fe oo i kK ok ok kK ok KK Kok KKK 


Beey,OGU WANT THE CONTROL OPTION? 


PieomeoheLUDES THE RICCATI SQUATION PROGRAM (RICATI) . 
eevot so ho COMPLET&E. DO YOJ WANT TO RUN LINCON AGAIN? 


PeieON IS NOW TERWINATED. 


mie weenmUr=sr Sesults £5" both options presented in 


(i) 


Appendix F indicates that a st32dy stat2= gain matrix has 


becr achizved. 


B) 





See DIisCRELE PINE KALMAN FILTER. (KALMAN) 

The discrete Kaiman filt2r gain matrix, G(k), is 
calculated by this program. The theoretical concepts 
involved in the davelopment sf the computer codes are 
eeceaped by Sege {[ 5]. A bris— deveiopment of the discrete 
Kalman filter is included her as an aid in the use cf 
KALMAN. Additional r2ference materials can be found in the 


Di bliograph y. 


Kalman fLiwer ing is 23 method of obtaining 
mMinimum-variance estimates aS Signais fron noisy 
measurements. The ijiscrete Kalman filter provides state 


estimates for the following systam 

x(k) = Sx(k-1) +4u(k-1) +Cw(k- 1) (4. D-1) 
with the discrete linzar observations 

Z(K) = Hx(k)+¥ (k) (4, D-2) 
where 


Bene Ne) State vestor at the tine “*(k) 


I> 
ul 


tO 
tt 


mie wae RONSINGULAL State trarsition metrix 


Memon cist iroence Eransiiizon of distributior matrix 


md 
i 





A = the 


the 


lx 
it 


the 


IN 
iT) 


Te) 
iT 
i+ 
DS 
o 


i< 
ut 
ct 
a 
Mm 


ic 
il 
ct 
1 
@ 


Mop Comerot diwceri bution matrix 

rx1 disturbance of syst2n randon input vector 
mx1 measurement vector 

mx1 measurement or observation vector 

mx1 measur2meant noise vector 

pel con<~rol or test signal vector 


discrete-time index (k = 0,1,...) 


The optimal filtered estimat>= of x(k), denoted 


XK (k|k), 
x(k |k-1) 
and 
Z(k{k) = 


Gor k = 


BK), is 
G (k) =P (k 
P(k[k) = 
P(k| k-1) 


where 


2s given by the recursive relations 


= bx (k-1|k-1) Yate 


E(k]k-1) + 3(k){z (k) -He (k] k-1) ] 


O, 1l,o.., Where x(0,0) = 0. The Kalman gain 


am nxm matrix which is specified is 

= 
|k-1H CHB (KI -1)d +R (kK) J-2 
[I-G(k) Hye (k] k-1) 


z 
= $p(k-1{k-1) 2 +9 (k,k-1) 


(i223) 


(4. D-4) 


iets 2 


(4. D=-5) 


(4. D-6) 


(4. D-7) 





I = the nxn identity natrix 


P(0|0) = P(C), the initial condition matrix 
x(k |k-1) = the single-stage prediction of x(k) 
x(k |k) = the filterei estimate of x(k) 


ls (w (k) wk) ME the nxn. 


Q= Sotaceance Nacrix of the 
Fandom inpat 
R = the mxn covarianc2 matrix of the measurement noise 
2 , ae ‘ 
E{w(k)w(k) }, the nean-square magnitude of the perturbation 


acceleration matrix 


1. Input Requireanents 

This projran comput3s the recurrence Eqs. 
(4.D-5), (4. D-5) ani (4.D-7) For. @ Speczized numbsce of 
meerat ions NP and prints 3(k) 25 2 Eunction of k. 

The required inputs are 

i} the seransition matrix (9) 

Zy the distributson matrix (P) 

3) E(w (k)w(k)" } 

4) the observation matrix (H) 

5) «he number of peints to be performed (NP) 


amc 


aeeitem = 2t gal Gond2-152 matrix [Pp (0{90) j. 


ui 
i 





Specific 
case. 
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~- 
— 
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by the user are 
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fx 
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Fauci feast Gletes 
RUD eleatied22 § OM 
eft Ff WGN weK 
COrmm mw IOS 
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WW Phd HaOmmamHOr 
dH OR FEU Ad 
HUH beet Py 


section 
Commands 
0.0] x(k) 
Oe 0 
of time poi 


Mert. 4 


Micon G 


thes 


= {1.0 


OHIGrHMOmOHOMRNO -» 


3 = fr, 
kal lon Se Se ca a 


MATRIX 


GAIN 


Wome ti (uP TO 8) . 


Oren 


Fk 
pw 


Nite, Simm i 


Ott AMABAHOZV 
=O = Oflchi<ct-o O 


VO MM MS MmsaHuW 


HHA AHOHREA A 


HONR AZ MEIMEmMOnM 


MU et gHO DUH MAatoO 
UAMMNVUFONAOAEG 
ft 


number 


Efw(k)w(k)"} = 4.0 
reo nh 


x (kK) 
R = (5.0] 


Zz (k) 
P (0) 


s 
= 


The following system is to b 


problem. 
The 
i 


iC_U 


rr 
ix 


B 


ame fn 
NOW, 


K 
Z 





ENTER T ELEMENTS OF THE TRANSITON MAIRIX--PHI 
PHI(1,1)= 
i 
PHI(1,2)= 
Bume(2, 1)= 
0.0 i) 
PHI(2,2)= 
THE PHI MATRIX (TRANSITISN MATRIX) 
WeOO r+ 0 0 5.00J0E-01 
Oe) 1.008 +00 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
nr 
ENTER THE DIMENSIIN OF THE 23ANDOM INPUT VECTOR. 
ENTER THE ELEMENTS OF THE DISTRIBUTION MATRIX--GAMMA. 
cases oT) = 


: GAMMA(2,1)= 


faeeonctA MATRIX (DISTRIBUTION MATRIX) 


eee 5-0 
5. OOE-0 1 


DOeYOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


Mitek he ELEMENTS OF THES YEAN=-SQUARE MAGNITUDS OF THE 
PERTURBATION ACCELERAILTION MAITRIX--W. 

W(1,1)= 

THE W MATRIX (eae lee eS PRONE? Or OF THE 

Pomeeon PATTON ACLE LERATION MATRIX) 
4. O00E+00 


Ome Welge TO CHANGE ANY ELEMENT OF THE MATRIX? 
Bvlee THE NUMBER IF OUTPUTS. 


meee > Hs BLEMENTS OF THE OBSERVATION MATRIX--d. 
H(1,1)= 


sete)! = 


mis HeafAT REX (OBS ERVATION 4AFRIX) 
1. 00E+00 0.0 
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DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


n 
ENTER THE ELEMENTS OF THE MEASUREMENT NOISE COVARIANCE 
MATRIX--R. 
a R(1,1)= 
THE R MATRIX (MEASUREMENT NOISE COVARIANCE MATRIX) 
5. 00E+00 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
al 
ENTER THE NUMBER JF THE POINTS TO BE PERFORMED. 
(<10 0) 
20 
ENTER THE ELEMENIS OF THE INITIAL CONDITION MATRIX--P. 
P(1,1) = 
1000¢ : 
P (1,2) = 
Py 2, 1) = 
B42e2) = 
1000. : 
THE P MATRIX (INITIAL CONDITION MATRIX) 
1. 00E+0 3 OeK0 
; 1,032 +03 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
al 


K= 0 

GAINS 

9. 958-01 0.9 
K= 1 

GAINS 

9.381E-01 1.922+00 
K= 2 

GAINS 

8. 295-01 9.98E-01 
K = 3 

GAINS 

7.03E-01 6.238-01 
K= 4 

GAINS 

6.132-01 4.52E-01 





Bley PAO 


Seo 0 | 


Baase =O | 


Sarco: - O01 


55 aye) 


See he 


Se 2 ve 


Cry] ye @ 


See eo 


Beeler 0 


Sac c Ud 


Bae a0 4 


™ 
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© 
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MG eDOO 


<7 





K = 19 
GAINS 
4, 86E-0 1 Sore ee 0 
K = 20 
GAINS 
4.86E-01 So 41a 07 
THIS CONCLUDES THE DISCRETE XALMAN FILTER PROGRAM 
(KALMAN). 


ANALYSTS IS COMPLETE. DO YOJ WANT TO RUN LINCON AGAIN? 
LINCON IS NOW TERMINATED. 


The computer results ar2 oresented in Appendix S$. 
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Ee. STATE VARIABLE FEEDBACK PROSRAM (STVAR) 


Given the linear time-invariant svsten 


Xx (t) 


Ax (t) + bu(*) 
a(t) = Kic(t) - k x(t) J 


o'x(t) 


y (t) 


the following can be performei by STVAR 


(4. E-1) 


(4. B-2) 


mppece Hellaermen 9E the plan: tianfer function, Y¥{s)/U(s) 


Wempoye Getifang a ficticiodus c¢ vector the internal 


BMS Hore rinct one can be calculated, X; (s)/U(s); for 


example, if X3 (s) 70 (5) wo decered theme Mates xis 


aopeeeeeca Wate cCs =) tate sdtl Othse ¢ Sebemhents squat 


emer te Gclez4s (5) ots d2Sirsd, calcuiars 


X, (s)/U (s) ani X,(s)/U(5) and divide he two 


[mec] | cile zon 2f the closed-loop transfer function, 


£(sS)/R8 (Ss) 


(iy caiculet=on of the Peedhacwes 2E21SiLe- -funct 20n, 
eq (S) 

[meson 2 Gesizei closed-id9p <transf2r Function, the 
Gon trailer Jase ana eam@back Coetiiczents zon 
media Opie, £6 Heq (S) eel be crtalculateda: the 





feedforward gain is salacted so that a zero steady 
State @rror cssulits fron a step input; ane designer 
who wishes other conditions must rescaie the gain 
and feedback coefficients appsopriacedy; ZOE 
exampie, if it is desired to hav? K = 1.0 but it is 
calculated as K = 2.0 with feedback coefficients of 


ho —sOm>, K = O.O0Mand k = 145, themprecedure “o 


2.0 
k = -=--- (0.5 0.0 1.5] 
1.0 


eeeleo 2Oe0 3404 


MiepeoL the iinfocnation necassa 


it 
KS 
rh 
O 
| 
ct 
om 
@m 
& 
in 
iD 
| 
ct 
O 
Nn 
O 
j~ 
< 
(iD 


state variable feedback probl2ms is presented ina this 
Section. However, th2 theoretical concepts involved in the 
development of the tomputer codes are fully described by 
Melsa [ 7]. 

Maem oasis =nput contains the problem identification, 


Metrices A and b and the order sf the plant, nn. These four 


Yl 


Mepues are requited reaqariless 9£ what opan- or closed-icop 
calculations are to b2 made. 

Poms oOent Siva R VerlLfiss =he controllability of <he 
meee. Th-ee Controllability conditions are possible 


Momecomoplece controllability 
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(2) numerically ancontrollable 


and 
fey) eunecnt ror ability 
Contollability arises when the controilability matrix 
E=[ b Ab Ath... AN-!b] | (4. 3-4) 
is nonsingular, Leo. jem cte ¥ 2. Even if the matrix is 


Demsoragular problems may aris® if it is difficuit to invert. 
To check this STVAR multiplies the matrix by its calculated 
inverse. The result should b2 the identity matrix. The 
actual matrix product is compar2d with tha identity watrix 
to provide a measure of uncontroliability. If the maxinun 
Merwe of deviation 1s not negligible, he plant is 
identified as numeric3 lly uncontrollable, A deviation larger 
Meemelo-s +O 10-5 h32s bean found to indicate difficulty by 
Melsa { 1]. 


yst2m is identified as b 


Ty) 


IC 
in 
iD 
Ue 
on 
D 
x 
yw 
| 
10) 
es 
tH 
mm 
ct 
iJ 
(ip 
OF 


ang 


DemeentrOliable, all dpen- ani closed-loso calculations are 


There are three possible closed-loop options available, 


iJ 
fas 


one for analysis purp Phewoche=> two for da si. 


O 


= 


Mm 


Sere ney oee 


Q 


n 
Pieposecs. After choosing the analysis option, and suppiyin 


PMyTR Wlth “he fteedforwerd gain K and ‘feedback cceificient 





matrix Ky the program calculates the closed-loop 
@aaractsristic polynomial ani the numerator of the 
equivalent feedback transfer function. 

The two design options are used t5 calculate the 
controller gain and the feedbask coefficients necessary for 
a desired Closed-Loop charazsteristic polynomial. The 
polynomial may be antered in either polynemial form or 
mactored forh. 

1. Terminal Session Examples 

HeeseesectiOn Contains a2 third order system for 
analysis by STVAR as present2d by Melsa = 1]. The state 


variable representation of the plant is given by 


e ~1.0 120 O20 O20 
mee 02002=~«~i ti(ié‘iCS Oe x CE) + 0.0 (4) 
eo ~=3.0 Ws - 0 
a I uL 4 
Weep) = f 1.0 1.0 CeO ei x=) 


{2s 
ct 
O 
th 
a) 
He) 
fru 
ct 
am 
(D 


Bee the Open-Loop cast, it is desire 
meeemtaai cranster function X (3s) /U{s) and the piént transfer 
Mec lon Y(s)/U(s). A fisticious matrix, necessary to find 


jmmemeneernal transtec function, is «hen 


2 = (70) 50) P nn Be | 





closed-lcop 
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values of the feedforward gain and the 
reguir2i 


coefficients 
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hla tole, 
NE mHO- 
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CART Xeak 


Cheeses ft el 


ELEMENTS OF THE PLANT 


HE 


Poeun LHS ORDER OF 


NOW, 
A (1,2) = 
A (1,3) 
A (2,1) 
pee » 2) 
A (2,3) 
A (3,1) 
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Tit A “MATRIX (PLANT asain 
Ore 0 0.9 te JUE+U0 
0.0 - 3. OOE+00 Jeo 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
n 
B(iys Lubec oe OF THE .CONTROL VECTIR--B. 
0 
B(2)= 
0 
B(3) = 
WE 
Lak B ieee (CONTROL MATRIX) 
0.0 
1. 0OOE+00 


DO YOU WANT TO CHANGE ANY ELEMENT OF TEE MATRIX? 


RE AK FE Ke KC A RK eae EE Ce RC oe ee a a he he EK eC 2 ee eC ae Ke ee eK Ce ee KOK 2K Xe 


OPEN-LCOP CALCULATIONS 


DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF S 
3. OOE+00 3.00E+00 1.)90E+00 1.00E+00 
THe ROOTS ARE REAL PART IMAG. PART 
e i -~1,73 F499 
0.0 - J3E+00 
-1.00E+00 0.0 
D@meroeu Have A FICTICTOUS OUTPUT VECTOR TO ENTER? 
y 
arOne Migeesoeee Ni > OF THE PECTIGCLIOUS OUTPUT VECTOR--C. 
0 
C(2)= 
0 
a5) = 
We (3) 
1g as pee CPs TTe1OUS Out 2UL VECTOR) 
OE 
1. COE+00 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
nh 


Meco OR COLE PEOLENTS = EN ASCENDING POWERS OF S 
0.0 ers 0 0 le) 0 +0 0 
56 





mas ROOGIS ARE R 
DO YOU HAVE ANOTHER FI 


THE ELEMENTS OF THE TRUE 


(OUTPUT VECTOR) 


© 


Sea PART 


0.0 
MATRIX [TO ENTER? 


SUPER serok-=C. 


DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


Mile nat OR COEFPFIOTIENES - 
. 00E+00 1.008+00 


THE ROOTS ARE REAL PART 
-2.,00E+00 
THIS PROGRAM IS CAPABLE OF 2ER 
OF CLOSED-LOOP CALCULATIONS. O 
MODE. THE OTHER IWO FOR DESIGN 
ENTER ONE OF THE FOLLOWINGS “OD 
1) -- FOR THE ANALYSIS 49 
2) P -- FOR THE DESIGN MODE 
CLOSED-LOOP CHARACT 
3) F -- FOR TYE DESIGN MIDE 
CLOSED-LOOP CHARACT 
GYOSED=LOOP CALCJ LATIONS 
KEY = Pp K 36K KK 
ENTER THE DESIRED CLOSED-L92P 
COEFFICIENTS IN ASCENDING PIWE 
YOUR HIGHEST ORDER COEFFICIENT 
POLY (1) = 
POLY (2) = 
POLY (3) = 
POLY (4) = 
THE NUMERATOR OF H-EQUIVALENT 
IN ASCENDING POWERS JF S| 
5. OOF-01 1.508+00 1.505 
THE ROOTS ARE REAL PART 
-5.008-01 
-5.00E-01 


e any 


tY Gj Ov) 


IN ASCENDING POWERS OF S 
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ERISTIC POLYNOMIAL 
2 A VALUE OF ONE. 





mis FE EPBACK COEFFICIENTS | 
B. OO8-01 0.) 1.30£+00 
THE GAIN = 2-000000E+00 
CLOSED-LOOP UCHARACTERISPIC POLYNOMIAL - 
IN ASCENDING POWERS OF 
&. OOE+00 6.00E+00 4.ODE+00 1.00E+00 
THE ROOTS ARE REAL PART IMAG. PART 
-2.00Et00 0.0 
-~1.00E+00 -1.00EF+00 
-1.00E+990 1.00E+00 
MAXIMUM NORMALIZED ERROR = Q.QO0 


DO YOU WANT TO RUN ANOTHER MODE IN STVAR? 


nN 
mio CUNCLUDES TAE STATE VARIABLE FEEDBACK PROGRAM 
(STVAR) . 
Pee ysis [5 COMPLETE. DO YOU WANT TO RUN LINCON AGAIN? 
n 


LINCON IS NOW TERMINATED. 


PRE COMmOutste wmesults £56 Ehis problem is presented 
in Appendix HH. The open-loop portion indicates that the 
PWecem tS controllable, since there is no indication of 
fimeontrollability, with thea desired internal transfer 


mercetion of 


Y(s S+2 


0 (Ss) S3t s2+35+3 


Mio wehoOccO-1 559) PObtion, Using the polynomial d2sign 


mode, shows shac the feedforward aqain and feedback 





coefficient matrix rayuired t9 yield the iesired ciosed-loop 


Srensfer function, Y(s)/R(s)}, are 


K = 2.0 


Or 
\O 





pees 





Eee LUENBERGER OBSERVER PROGRAM (LUEN) 

Wet a “pacticular slosed-loop transfar function is 
desired and some of the states ar? inaccessible, LUEN can be 
used to design a Luanbergyer Jbserver. Poa example, 22 gq 
measurements are state variables, an observer of reduced 
dimensions can be designed ¢9 32stimate only those states 
which are not measur3sd. The state estimates generated by an 
observer can be used as PEMDUEt = RESENatioOn ae, a 
eenctroller [{ 8}. The block diagram presented in Fig. 4-1 
represents the general form of the system when a compensator 


is placed in the f2edback path. 







r(t).Jcontroller 





observer 


Fig. 4.1. Luenberger Observer Block Diagram 
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The plant is characterized oy the stat2 and mnéeasursment 


equations 
X(t) = Ax(t) + Bu(t (4. F-1) 
z(t) = Cx (t) (4. F-2) 
z*(t) =o x(t) (4. P= 3) 


An observer can b2 designed that generates an estimate y 
which converges to the state x as time becomes large — 9]. 


A linear controller is design2d as 
u(t) = K(xr(t) - k' x(t) ] (4. P=) 


2-n which all states can be ma3asured. Replacing the true 


State with the estimate yields 


u(t) = KE) - kT y(t) J (4. F-5) 
where 
Rye = b(t) + g z ft (4. F-6) 


; T ; | 
AS time increases k y(t) will aporoach k x(t). Since the 

reduced-ordeéer observer 2Stinates only the inaccesible 

States, the control law contains measured states waers 


[eee ble and ObServect Estimates for the other states. 
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and 


Eoumerse in using LUEN, the £>llowing are defined: 


X(*) 


(> 


io 


IA 


i> 


5° 


y ( 


he 


) 


HN 


i 


n-element column st2zte vector 


plant input 


q-vector of systen measurements 


output variable t) be controlled 


plant M2trix (nm x ny) 


distribution matrix (q x m) 


[assumptions: (1) J $n and (2) € has rank qj 


output natrix 


fOECERG: Lane ca on 


herimouw ard OG seGa2cOLlor gain 


EeeapacKk CoOstitessneemmatrux{n x m) 


@pecervec feedback cosificient matrix 


(iD 
er, 
Le] 
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Q 
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From theory, the observer is described by 

y(t) = EFy(t) + 3,2(%) + G,u(t) (ic ie 3) 
where 

F = observer eigenvalu2 matrix 
and 

G, and G, = observ2r gain matrices. 


Z'(s)/R(s) can be solved by the following procedures 

(1) select Z'(s)/7R (s) and solve for K and k by using the 
STVAR progran 

(2) use the OBS CON projram mie) calculaze2 the 
observability index the observer is designed tod 
have a dimension greater than or equal to the 
observability ind=x minus one. 

fey select the F eigenvalu3s; these should not equal 
emocenwon A, JEeVEoDuUSLY S2lculated by SIVAR. 


(4) use LUEN to calculate G,, G,, h and gq. 


The 2xample ocesented h2r2 is taken from Desjardins 


mee l47 ( 2]. Commands entered by the user are in Lower 


rie 





Gese. However, due to the spasial natures of the example, 


comments have been added tds the terminal session for 


Silarity. Given the fourth ordec plant 


COmn..0. GOmoumOL 0" ramon 
<, 020 9:0 1.9 0:0 a. 4 0:0 se 
xu To = x u 
x(™) VO Wee wae ee 0.0 
O15. (253,20 1:0 
| Jj i y | 
with 
z*(t) = [10. 20. 0.0 0.0] x(t) 


Because x, and x, are the only neasureable states 
iS SO OMe Nec ren) SoG) ey 
= Xx 
zt) 2.0 1 tao 020 ~ 


The closed-loop transfer function to be achieved is shosen 


to be 
Z* (s) _ 1 
R(s) $44+55341752+285+20 
lincon 
reese UTLON BEGINS... 
meeon CONSISTS OF [THE FOLLOWING SUBPROGRAMS: 
Eeeoec MATRIX MANIPULATION - <BASUMAT> 
MmeeroyvAl TIME RESPONSE - CRIRESP> 
STATE VARIABLE FEEDBACK =- <5IVARD 
Se PenOLIASITLEITY AND OBSERVABILITY - <OBSCOND>D 
Met oeRGER OBSERVER - <LUEN> 
Seem al CONTROL/SS&LMAN FILTERS ~- <RICAIiI>D 
Ptesenel & EEMh CAALMAN FILTER = <KALMAN>D 
CPTTHAL CONTROL - <OPICON> 
Pee eaL FRACTION EXPANSION - <PREFEXP> 
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THE NAME BETWEEN THE 


Pei an THE GRDER JF THE SYSTEM (UP TO 8). 


NOW, 


uy 


Pee DLHE ELEMENTS OF THE PLANT MATRIX-<-A 


A(1,1)= 


A(1, 2) = 


A (1,3) = 


A(1,4) = 


A(2, 1) = 


A (2,2) = 


A (2, 3) = 


A (2,4) = 


A (3, 1) = 


Be 5, 2) = 


A (3, 3) = 


A (3,4) = 


Ree) = 


A(4, 2) = 


~ @ 


Rea, 3) = 


Je 


A (4,4) = 


ue 





THE A MATRIX (PLANT MATRIX) 


0.0 1,.00E+00 0.) 0.0 
0:0 0.0 1,00E+00 06.9 
20 0.) Hod. 1, 00E+00 
0.0 -1.50E+01 -2.30E+01 -9.09E+00 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 
n 
ENTER THE ELEMENTS OF THE CONTROL VECTOR--B. 
B(1)= 
0 (1) 
B(2)= 
0 ( 2) 
B (3) = 
0 ( 3) 
B(4) = 
a (4) 
THE B MATRIX (CINTROL MATRIX) 
0.0 
0:0 


0.0 
1. QOE+00 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


ME KK RI Ke HK EE EK RK Re HK Ke CK He eK Re ie eK eK ie A OK eK EK KK kK OK OK 


OrPEN>LOOP CALCULATIONS 
me onl wat OR COEPY ECLTENTS - IN ASCENDING POWERS OF 5S 


O30 1.50E+01 2.302401 Y9.O0E+00 1.00E+)90 
THE ROOTS ARE REAL PART IMAG. PART 
=3,00Et00 0-0 
-5.00E+00 0.0 
=Peooeeoo 0.0 
0.0 On 
DO YOU HAVE A FIZSTIICIOUS OUTPUT VECTOR [0 ENTER? 
n 
Breileem THE ELEMENTS OF THE TRUE OUTSUT VECTOR--C. 
Ci1i)= 
25). a 
C(2)= 
reo! 
C(3)= 
0 (3) 
C(4)= 
Giceee barn tx (OUTPUT VECTOR) 
2. 00FE+01 
1. 0OE+01 
O50 
0.0 
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HOmroOu ONAN. TO CHANGE ANY ELEMENT OCF FHE MATRIX? 


n 
NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 
2- COE+01 1.008401 
THE ROOTS ARE REAL PART IMAG. PART 
-2.00E+00 0.0 
THIS PROGRAM IS TCAPABLE OF PERFORMING THREE TYPES 
OF CLOSED-LOOP CALCULATIONS. ONE TYPt FOR THE ANALYSIS 
MODE. THE OTHER [dO FOR DESIGN. 
ENTER ONE OF THE FOLLOWINS YODES: 
ieee -—- BOR STHE ANALYSIS MODE 
Zee -~- Pores DESTGNSUODS WITHVTHE UNEFACTORED 
CLOSED-LOOP CHARACIERISTIC POLYNOMIAL 
3) F -=- FOR THE DESIGN WIDE WITH THE FACTORED 
CLOSED-LOOP CHARACTERISTIC POLYNOMIAL 
p 
SLOSSD-LOOP CALC I LATIONS 
KEY = P fe Yee oe oe 
ENTER THE DESIRED CLOSED-LOOP CHARACTERISTIC POLYNOMIAL 
COEFFICIENTS IN ASCENDING POWERS OF S. 
YOUR HIGHEST ORDER COEFFICIENT MUST BE A VALUE OF ONE 
POLY (1) = 
POLY (2) = 
POLY (3) = 
a7. 
POLY (4) = 
POLY (5) = 
THE NUMERATOR OF See = 
IN ASCENDING POWERS ESS 
2e QCOE+0O1 1.3 0E+01 -5.90E+00 -3.00F£+00 
THE ROOTS ARE REAL PART IMAG. OART 
-1.0%E+00 0.0 
ec PR a0 0 A 8) 


Poioet seUBACK COEFFICIENTS 
peeve eon 1.508401 =—6.008 400 =—3.0058+00 
THE GAIN 1.00000 0E+09 


Pie owOseb-LOOP CHARACTERISTIC POLYNOMIAL - 
MiP asSeC SNDEING POWERS JF S 


Peeocetoie 2.508*0t T.70E401 6.00E+00 1.00E+00 


Uz 


V7 





THE ROOTS ARE REAL PART IMAG. PART 
-1.008+90 -2.00E+)9 
-1.00E+00 2.00E+00 
-2.-00E+00 0.0 
-2.0U0E+00 0.0 
MAXIMUM NORMALIZED ERROR = 0O.0 


DO YOU WANT TO RUN ANOTHER WIDE IN STVAR? 


THIS CONCLUDES THE STATE VARIABLE FEEDBACK PROGRAM 

(STVAR) . 
COMMENT: jrcasults indicate that the system is 
GommrOllebie, that th: ee SEA wee be = 3.0, 
aoe, cisco and 0.0, that the values ofek are 20.0, 
eres Oe0 2290 —-=3.0 and that K is equak.to unity; 
the observability ind2x is next needed to design +hea 
observer 


ANALYSIS IS COMPLETE. DO YOJ WANT TO RUN LINCON AGAIN? 


Y 

LINCON CONSISTS OF [HE FOLLOWINS SUBPROGRAMS: 
BASIC MATRIX MANIPULAIION - <BASMATD>D 
RATIONAL TIME RESPONSE =- CRIRESP> 
STATE VARIABLE FEEDBACK - <5SIVAR> 
GeurROeLIASTLITY AND OBSERVABLILLITY ~ <OBSCON>D 
LUENBERGER OBSERVER - <LUEND 
OPTIMAL CONTROL/KALMAN FILIERS - <RICATIID> 
DISCRETE TIME SALMAN FILTER - <KALMAND 
OPTIMAL CONTROL - <OPICON>D 
PARTIAL FRACTION EXPANSION - <CPRFEXPD 
ROOTS OF A POLYNOMIAL - <RODIS> 

TO USE ONE Of THE SUBPROGRAMS ENTER THE NAME BETWEEN THE 

SYMeOLS < >. 

obscon 
OBSCON DETERMINES THE OBSERVABILITY INDEX AND 
@GOnNT ROL LABILITY DF A SYSTEM. 
mibeot ENTER THE PROBLEM IDENTIFICATION 
MenOoT TO EXCEED 2) CHRRASEERS*). 

thesis example 

; Meng) Sluesn THE ORDER QF THE SYSTEM (UP TO 8}. 
elon LHe SLEMENSTS OF THE PLANT MATRIX--A. 

0 A(1,1)= 

1. A (1,2) = 

9 A(1,3)= 

0 A(1,4)= 

5 A(2,1)= 
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(2; 2) 2 


0) 

1. A(2,3)= 

0 A (2,4) = 

5 A(3,1)= 

» 402) 

0 A(3,3)= 

1. A(3,4)= 
A(4, 1)= 

0 (4, 1) 

A(4,2) = 

a 
A(4,3)= 

oe , 3) 

_o Seti 
aie A mee l RIK (PLANT WATRIX) 

Q.0 1.902409 oae0 ae) 

O70 Q.9 Vs 0C0l #005) =o 

ANG, 0.0 " 2 UO 00 

O70 - 1,.50EF+01 -2.30E+01 -9.00E+00 
meer, OU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 

1g! 

1 Bese Lae NUMBER JF COLUMNS OF THE 38 MAIRIX. 
Pei ores Oe SLEMENTS OF THE B MATRIX. 
B(1,1)= 

0 

0 Ex(27etj-= 

5 B(3,1)= 

i. B(4,1)= | 
foe o MAPRIX 

6.0 

0.0 

ome) 

1. 00E+00 


Some OUmewewe 10 CHANGE ANY BLENENT OF THES MATRIX? 
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ENTER THE NUMBER JF OUTPUTS. 


MATRIX. 


(S 
‘a 


THE ELEMENTS OF THE 


R 
7 


C (1,2) = 


C(1,3)= 


C(1,4) = 


C(2,1)= 


C(2,2) = 


C(2,3)= 


C (2,4) 


0 
Q 


J 
QO. 
DO YOU WANT TO CHANGE ANY ELEMENT OF 


THE C MATRIX 
OE+O00 0.0 
1. 00E+00 


0 
0 


1 
0 


Doane tee Ke 


3 


OBSERV ABILITY INDEX 
ME Ol ol el 
AE OK eK OK KK KK OK eK Kk Xe 


2 HE KE CK Ne ae Oe oe OK EK OK OK 


LE 
mK XK 
? 


Mt = 


ae 


eo 


DO YOU WANT TO RJ 


NDex AND 


4 


erp ete 


a) aya 
ae oul 
mm ra) 
OW 


“rhc tn 


m+) a) a) 
a aw 


«tl «1O 


ON at.G 
WNP >t 
eA 
Ae 
Matha 
W Wert) 
Mma) 
ta tin 

cs 0) 

MNeeDD 
eee ft HHO 
Ec a @ ws 
Zhe> > 
eo & 
=) W1Q 
=f Wom 
Ole. Fieri 
UM o~w” 


DO YOU WANT TO RUN LINCON AGAIN? 


omc ee be 1.5. 


ANALYSIS 


U 

Ss 

iE 

R 

fio O53 CON? 

= sneeGar i> 
<KALMAN>D 

Rep yeeaee 2 

R THE 


MAR 


B 
E 
V 
L 
= 


(HV MmFHH Mm 
ar FAInmMaA fr 
Hl mV tao oO 
Ss VV > WAMA = aM 
Oz (MDOHHZOVES 
HO I Aart JOW eG 
HH tMnvyvY HOW 1G 
Orn MM mhurie G3) 
fueQTWIO’rD Aledo 
= 22 O Auer 
WMOAA MAY MEIN 
I$ Qu Mu f23 & Fa ert fst 
F444 U9 Fr ecrit> bi 3 orm 
Mii MOND wmaeWw 
Pfaetm HMO AO Ms 
Oe 


VSS OAOHMHOUME 
HHH CoH Mbit 

NM HHMOMGO HO im thy 
He FIO Of Oe 
UY ett} C) fz) yn 
Sati IGA AHAOR 
O = Ofidthidce ZV 
OVO MMB mBHNO 
HH HESHOHHE W 
ZN Hea SHEN E GOA 
Ot tHODUHMAtONO 
OMMmNUAOAOMAmDM”A 
: 5 


HH 
mm 


Lluen 


3) 
EUW) 
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ENTER THE 
ENTER THE 
ENTER THE 


ENTER THE 


Be 2) = 
A(2,3)= 
A(2,4)= 
a3, 1) = 
A(3,2)= 
Rt: 3) = 
A(3,4)= 
A(4,1)= 
A(4,2)= 


ths) 


A (4,3) = 


A(u,4)= 


Ntd Fi WM 


ED T) DESIGN LUENBERGER OBSERVERS TO 

GIVEN ZLOSED-LO3P TRANSFER FUNCTION WHEN SOME 
ABLES ARE INACCESSIBLE. 

R THE PROBLEM IDENTIFICATION 

CEED 29 CHARACTERS). 

e 


ORDERS OE THE SYSTEM (UP TO &). 
NUMBER JF MEASUREMENTS (UP TO 8). 
ORDER OF THE OBSERVER (UP TO 8). 


BEOWENTS OF THE PLANT" MATRIX--A. 
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20. 


THE A MATRIX (PLANT MATRIX) 
0.0 PaOureu0" J.) On 
0.0 0.9 1.J0E+00 0.0 
0. 0.9 J.J V0 E00 
0.0 SY) SUG S OSS Aa EOC SO) 1) ote) 10)19) tie 210118. 
DO YOU WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


ton hoes LeWeNtS OF THE DISTRIBUTION MATRIX--B. 


B(1)= 

B(2)= 

B(3)= 

B(4) = 

THE ; MATRIX (DISTRIBUTION MATRIX) 
0.0 
0:0 


Q00E +00 
PO YOU AANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


Seek SHE ELEMENTS OF THE OUTPUT MATRIX--C. 
Cc (1,1) 


C(1,2)= 
C1, 3) = 
C(1,4)= 
C(2,1)= 
C(2,2)= 


C(2,3)= 


Oe MAT ALA: 


ENTER THE DESIRED FEEDBACK COEFFICIENTS. 
PDBK COEFF (1)= 


2) 





hook COEFF (2)= 


13. 


FDBK COEFF (3)= 


= © 


FDBK COEFF (4)= 


Ss 


Poop soln kD FEEDBACK COEFFICIENTS 


O CHANGE ANY ELEMENT OF THs MATRIX? 


CTroor 


FORM 


CLIENT 


COBL ED 


£ 
an 


ENTER THE REAL PART OF THE ROOT. 


> 3 Ape, 


MAGNITOIDE OF THE IMAGINARY JOT. 


ENTER T HE 


THE ROOT. 


ONTER THE REAL PART OF 


MAGNITODE OF THE IMAGINARY ROOT. 


Tt HE 


R 


FIM = 
me) 


OBSERVER EIGENVALUES 


KOR AK ERC TK ok He Ke Oe OK OK AK KK KK OK KK OK KK OK OK KK KK eK ok KKK KKK 


DavokR SLSENVALUE MATRIX) 
0E+00 


Be SAIN MATRIX) 


MATRIX) 


GAIN 


(OB SERV ER 


On ies: 


er 


Soa ica SOs ret 


SATOR 
COE +00 
bioese Pais LUENSERSER OBS 


a 
2 


OM Pe 
; 


c 


R DESIGN PROGRAM. 


yn 


Ne 


ol 
- 


aS 





DO YOU WANT TO RUN THE PROGRAM AGAIN? 


ANALYSIS IS COMPLETE. DO YOJ WANT TO RUN LINCON AGAIN? 


LINCON IS NOW TERMINATED. 


The results of STVAR, I2BSCON and LUEN are shown in 


Appendix I. From these th2 the observer is given as 


Meeoet0 Y(t) 9 8525 29.25 x, (t) | -3.9 
+ 


y(*) = + u(t) 
714.0 0.0) yy (t) | 929 060 | x, (td, as ey 
and 
c rs a | + a | 
a(t) = 1.0 ¢e(t) - (20.0 8.5) “8? ~ 61.0 0.0) 23°") 5 
4 Xo (oie ae Cell) 
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G. OPTIMAL CONTROL PROGRAM (OPTCON) 


Given the linear, time-invariant system represented as 
x(t) = Ax(t) ¢ Bu (t) (4. G-1) 


OPTCON will minimize the cost function 


N-| 
3(N) =1/2x (N)Qx (N) +172 4 [x(k) 2x (kK) + RU2 (k) ] (4. G-2) 
Rz0 


where 


State vector 


Ib< 
ul 


Q = measurement nois? covariance matrix (n x n) 
N = number cf tine intervals over which the sum is made 


R = random input (a scalar) 


A = plant matrix (nx n) 
Pee d2stribution matrix (n x 1) 
ane 
icy) = cCentrol (a scalar). 


The output of th2 program is the feedback gain natri 
Pech, When muis=ipli2i by the state vector, yields a scalar 
Seen. =Ol. The following recursive equations were derived 
Meng Gynamic programming, Starting at the terminal time and 


Me -kKing backwards. 


2 
P(k) = Y mre cx-1)Yer + Q + RA(K)A (Kk), B(0) = 9 (4. G-3) 
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Woy = } + Aa’ 9) VAT = 0 (4.G-4) 


ts 7] 

a’ (k) = -[AB(k- 25/74 Pk-1)4 + RJ, a (0) = 0 (4. 6-4) 

For Simplicity in programmniag, the following terms are 
defined: 


T 
terminal = 1/2x (N)Qx(N) 


Nel 
trajectcry = 172) x(k) Ox (<) 
det R:0 
fuel = 1/2 Rue (k) 

Rzo 


1. Teroinal Session Exampl:; 


Given the system and parameters described below find 


the discrete steady state gains for a sampls of 0.1. 


x (t) 0.0 a a | (2) 
x“ = a OT Gat 
= Tae) 80.0. > 1.) 
po Jj au 4 
Meio) | (1.0 
‘1,0 | Bien 
J 
Re= 1.0 


iomeOods= TON, cun the program fort 2 time interval of 
GQ. 


Pe acon 
Meee UTrtON BEGINS... 
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a 
HH 
a 
fxQ 
fxy 
= 
an 
AN fQ 
= aa) 
ee © 
ep) UO A fx 
= Ww KH = 
<c mM Fa ax 
ee OO aad a 
C2) ¥ USA A 
OA Hie Ay FQ 
OG EH ' as MM 
Auer VH fF a 
MSAAM Gf 
MWMMmA i'M wGw 
Nt W) iH Vv AAR 
MNhitfyH WwW VWwR 
(OV MF at on 
= RtnoMmaA fr (Ow) 


MlmVetz Ho O 
S$ VY PMIhlAS MUN 
Om tMDHHZOVS 
HO’ MAmMWAOH « 
FH wMNY HUW 
OmmRIOM whhHae © 
MaMBO et Octo 
= SION a 
MDOAAM~ALHV «MHA 
THUG leat st Ist 
Fs latte kid tt 
Mh MmNnad BMmy 
fudio; Sef eQom 
OS RMHNO OHH 
DAHWMmMMmAOm 
USS OAOR SHU Ae 
AHHH wHmet 
NG HMOGmO HOmMcdth, 
ies MGenthor)y O)fsi O 
Creo fr 
BES ait AG IQHAWOn 


Bs A 


QO = OhWatmade zAV 


VVO NGA asHnNnoO 
HHA AZHOHHE 


A O 
oa & 


WY) 
2WNbrsea2MHnHmOrA 


MINIMIZES THE THE FOLLOWING COSI FUNCTION 


J (N) =MIN (SUM (X(N) THQ¥X(N) ¢UT (N-1) *R*0 (N-1))) 


OPTCON 


it 
pa 
ae EA © 
~<O =~ — 
HHF4 Cheat 1H 
moO mm wh) 
Bm HM Ay 
t> WS moO 
= [mti + it 
fx] HH He 
mae A fr2 
Hea FIN A~ 
Se to + 
ON HR «2 
mH = 
<q ws i) 
Or am 
Kir QoQ 
mM MA A 
Am ma + 
moO mM 6 
fr] mm 
fH E+ fx 
mM Wh a) 
he 2H! 
Tat Oe N 
Fay HI tH 
HAO Mm 
Tal | ro, Me 
Hy DH + 
co 0 @ | 
=) Fma =— 
«x . a 4 
Mane ts 
OROPrPNnN ~~ 
Om maH Ay 
Mee = + 
me 2ma0 8 


PSI(K) =PHI+DEL*AT (K) 


(2) 


P(Q)= 


(3) P(K) =PSIT(K) *P(K-1) PST (K) +Q#R*A(K) #AT (K) 


thesis example 


CATION 


eee 
*). 


OVER WHICH 


IME INTERVALS (NSTAGE) 
N 
H 


MSE aD TO 8). 
T *%+# 


SUR NUCBER «ek 


MAT *** 


OG D4 


Q(1,2)5 


Q (2,1) 


O0(Z,2)= 
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fx] 

ase 
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= 
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MsOOH 
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Fie eD 
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AT 


ee Dae SAMPLE ENTERVAL--Df. 
ae ONCE R DL TN F-FORMAT *** 


tJ 
pee 
id 


NE OF THE FOLLOWING OPTIONS: 
RO. 
E 
Mi 


ae 
ae 
rH 
OQ 


MINAL+TRAJECTORY+FUEL, 
NAL+TRAJECTORY+ 0 


im oO 
tH 
WwW 
N 
ty 
es 


THE NUMBER 1 IF 
E, COST=PERMINAL+ ) ¢ FUEL, 


R 
[ 
O, CIST=TERMINAL + 
A E A AND 


gr" Q ar" OH 


R + 
IF YOU WANT TO RE TRI 
BUT NOT THE S, fR AO 
E 


i Tiha be os, 


rT 


0 
BMA 
ATRICE 
HE PHI AN 
CES, ENTE 
MATRIX--A. 


0 
CEs 
ENT 
DD 
R A 
A 


M 
ic 
iL 
iy 


"Wr 
Ha 0 
try 
ty 
tr 
td 
He 4 
ty 
= 
a 
O 
rij 


1.00B+00 0 
DO YOU WANT TO CHANGE ANY BLEMENT OF THE MATRIX? 


Mibeseeeones Of THE DISTRISUTION MATRIA--8, 


EB MATRIX (DISTRIBUTION MATRIX) 


DO OYO@U WANT TO CHANGE ANY ELEMENT OF THE MATRIX? 


THE PHI MATRIX ; 
1.00 5003#+00 1.0016572-01 
1.00 1667E-01 1.0050032+00 
THE DEL MATRIX 
5.00416 48-03 
1.00 16672-01 
MINIMIZATION OVER ALL STAGES 
g( 1 = AZ (J) *K(I) 
AT (1) = -2.462-02 
AT (2) = -1.98=-)1 





a2 





e 
EACI ON 
egy 1 


" 

noi 
= 
~~~ 
NO-N 

be inaatl 
~HE 
Detet 


20, 


wee mee eee eee ee ee Ne Ne ee” 





% fy 


that 


tyascats 
and -2.4. 


JaRAM (OPTCON) . 
ire. 


WANT TO RUN LINCON ASAIN? 
BL 48 


Appendix J, 
about 


de 


= 47 
91 


a 
— 


DO YOIJ 
shown in 
achiev 


af 2 


20S 


LINCON IS NOW TERMINATED. 
mae cesuits, 


AHALYSIS 15 COMPLETE. 


steady state ga 





Ve. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Although LINCON was acicten primarily as 2 
teaching/learning tool, it can still be guite useful to the 
practicing engineer for design and analysis problems. [t was 
written in modular form so that it could b2 easily nodi fied 
by the addition of subroutines. 

LINCON has been e¢xtensiv3ly tested in an advanced 
optimal estimation <cours2. Th interactive aspects proved 
highiy successful. Hopefully all the "bugs" have been 
eliminated. 

Memormeated Carlier, althougi the original intent of this 
thesis was to adapt )Jesjardins' version of Melsa's LINCON by 
femeng it interactive, LINCIN began t2 grow as other 
routines were added and/or extsisively nodified. 

OPTCON, THew Pel Mal Geemerohe program usindi==necur sive 
equetions derived fron dynamic ocogramming, is a new member 
momeene LIVCON family. Desjardicst KALMAN, the discret> time 
Kalman filter progr2an, underw2int considsrable progrimming 


changes before achieving its present forn. 
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Be RECOMMENDATIONS 

(1) The source program, LINCON FORTRAN Al, must now be 
passed from user to user and than compiied before it can be 
used. As can be seen in the terminal session examples, the 
program was invoked by typing "iincon". Aeculeliyv, this 2s 
an executive progran, LINCON EXEC A1, comprising of the 


following statements: 
Meee per eOSePRINTER (RECFM FA LRECL 133 BLOCK 133 
LOAD LINCON (START 


Meeweirst statément iefines 09 as the printer and permits it 
eemepoiny out 133 characters per line. Ths second statement 
invokes the compiled version of LINCON. It 1s recommended 
that LINCON be placed on 2 utility disk so that users may 
Men tO 2t instead of the current procedure. 

(2) Peweeimes Lt Can be exceedingly ditfticuit to 
Mmeper pret the tabular output of some of tha pregrams. [2 is 
recommended that a@ graphics pactage be developed for ARIRESP, 
RICATI, KALMAN a@nd OPICON. The package should be interactive 
Weenie tne Output being first displayed on the terminal scrze¢en 
mmmeenon allowing <=h2 user £5 choose the type of output, 


$.€., VERSATEC, TEKTRINIX cr print-pict. 
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(3) In general, the programs are limited to eighth-order 
problems. If the ne3d should arise to solve higher order 
systems, this linizsation may 52 removed by extending the 
appropriate dimension statements. The user must remember to 
alter the format statament pertaining to the output, either 
decreasing the significant figures or adding a "wrap-around" 
feature +0 overcom2 the printer linitations =e 33 


characters per lins. 


j-> 


(4) AS it is written, LUEN can only solve for a 
treduced-order observer. It is recommendel that the program 
be modified so that the user has the option of selecting a 


reduced-order observ2ar or an idsnotity observer. 


(ayeeeene ficticidus and real ec matrices of SIVAR are 


$ 


required +o kave thea dimensions 1 xn. It is recommended 


that tne program be nodified to 2ccept a dimension size of n 
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KPPENDEX 5 


PARTIAL FRACTION EXPANSION 


PARODO EN TDENETTETCATION: Pino LS 
eK RR KK HK KEK RK RR KK KK ARK KK KK KK KX 


aise feel L Sea sAIN 


EXAMPLE 
HK KKK KK RK KK KKK KKK 


7.0 C00 OO0E+ 
He He He He He KK He ke ee Ke eK eK Ke Ke KK eK KK ee KK KKK KK KK KKK KR KKK ae XK HK XK 
NUMERATOR COEFFICIENTS - IN ASCENDINS POWERS OF 
2.0000E +00 2.39 0OF+ 00 1. 0000E+00 100008400 
NUMERATOR ROOTS ARE 
REAL PART IMAG. PART 
0.0 -1.474213E+00 
0:0 1.414213E+00 
- 1.0009 00F+00 0.0 
Me Ke OK te ee He ee Ke eK Ke aK ke KK ke KK KK KK KK KK KK KK KK KKK KK KK KKK KKK 
THE INPUT FUNCTION 3AIN 
7.0000 00E+00 


He ee eK tk KK KK KKK RK KKK KKK KKK RK KKK KKK KKK K KKK KKK 
DENOMINATOR COEPFICIENTS - IN ASCENDING POWERS OF 5S 


22 000E+ 00 1.000 E+00 0.9 3.000E+00 1.Q000E +00 
DENOMINATOR RODTS ARE 
REAL PART IMAG. PART MUDGEE Eeie ye 
=3. 0 18559E+00 0.0 1 
Soe OOF eo I 0.0 1 
4.6580 98E-91 woes oe 7 0 1 
Pogo 0 98R=) J Se sO G42) Bm 1 
kk doko ak ok dok kok doo kok kkk ok ek KK oR RK a RR KK KKK KEK 
REseDiUE OATRIXe=- REAL PART 
9 8 108 7 2ZE+00 
Bee 37 20 | 
4. 36023 1E-01 
aeesiou. 3 12-01 
Pees MATRIX - IMAG. PART 
wate 
2.0614 06E+90 
-2.0614 06c+00 
ek tk dakota ok kk ec ko oi oko ook oe ak deo ek kek ke ok kok kkk kk Kk Kok 
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OPTIMAL CONTROL/CONTINJOUS KALNAN FILTER PROGRAM 
PROBLEM IDENTIFICATION: THESIS EXAMPLE 
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APPENDIX G 


DISCRETE TIME KALMAN FILTER PROGRAM 
PROBLEM IDENTIFICATION: Rieols Ex awe le 
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